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Abstract 

In this work we show that the latest LHC data on multiplicity moments 
C2 — C5 are well described by a two-step model in the form of a convolution 
of the Poisson distribution with energy-dependent source function. For the 
source function we take V Negative Binomial Distribution. No unexpected 
behavior of Negative Binomial Distribution parameter k is found. We give 
also predictions for the higher energies of 10 and 14 TeV. 



One of the widely discussed, yet unsolved, problems in high energy hadron 
scattering is the production mechanism of low and medium pp hadrons. Here 
perturbative QCD cannot be applied and one resorts to phenonemonologi- 
cal models and/or Monte Carlo generators [Tj. With the advent of any new 
hadron accelerator the quantities first studied are charged particle multiplic- 
ities. Indeed, the first physics LHC paper, the one by Alice collaboration [2] , 
dealt with the average multiplicity. By now data on multiplicity moments 
were published by Alice [5] and CMS [3|. It is therefore important to find 
if their behavior can be explained in terms of some simple phenomenologi- 
cal models or the widely used MC generators. The latter has been recently 
addressed in Refs.[5]. 

Throughout this paper we explore the well known observation, recently 
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recalled in Ref.[B], that multiparticle production can be described by the 
probability distribution P(n) which is a superposition of some unknown dis- 
tribution of sources F, and the Poisson distribution describing particle emis- 
sion from one source. This is a typical situation in many microscopic models 
of multiparticle production. For example in Dual Parton Model (DPM) par- 
ticles are emitted by chains spanned between the colliding protons (for review 
see Ref. 0). In that case there are many sources whose number increases with 
energy. Similarly in Quark-Gluon String Model (QGSM) emission proceeds 
from n cut-pomerons (Ref.[8] and references therein), and again is assumed 
to be Poissonian. 

Here we refrain from formulating a microscopic multiparticle production 
model and assume a simple phenomenological formula that captures, how- 
ever, the above mentioned physics of independent emissions encoded in DPM 
or QGSM: 



P{n) = I dtF{t)e~ mi ^. (1) 
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Here t is a fraction of the average multiplicity, and F(t) the distribution 
of sources that contribute fraction t to the multiplicity probability P(n). 
Normalization conditions require 

oo oo 

JdtF{t) = JdttF{t) = l. (2) 



There are two useful properties of Eq.Q that will be of importance 
throughout this paper. The first one is the fact that factorial moments of 
multiplicity distribution measure directly the F m+ i moments of the source: 



(77(77 - 1) (n - 2) ... (77 - m)) = n m+l J dtt m+1 F(t) = n m+l F m+1 . (3) 



Because of (J2J) average multiplicity 

(77) = 77. (4) 
Factorial moments can be expressed through scaled regular moments 

C (5) 

(77) 
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that have been measured at the LHC [21 H]- For the first five moments we 
have: 




(6) 



The second property of Eq.([T]) is that for large multiplicities (i.e. for 
large energies) it implies an approximate KNO (Koba, Nielsen and Olesen) 
scaling [2]. Indeed, in the limit n — > oo and fixed n/n one can apply the 
saddle point approximation to calculate dt integral in (fTl) leading to [6J 



KNO scaling says that function tp depends only on r = n/n. 

KNO scaling is seen approximately in the multiplicity distributions mea- 
sured at SPS and higher energies (see for review [10]) including the LHC 
[3J H]. This fact gives strong justification for formula ([TJ which also allows 
to give definite predictions for the violation of the KNO scaling. Originally 
KNO scaling has been derived assuming Feynman scaling [TTj . which states 
that the central rapidity density saturates at asymptotic energies. The latter 
is clearly not seen in the data (e.g. [2]), on the contrary central rapidity 
density grows as a power of energy. For example for \q\ < 0.5 |12j : 



where W = \/s. For constant n all multiplicity moments would be constant 
as well. 

Here we see the advantage of the convolution model ([T]) since it implies 
approximate KNO scaling also for energy dependent n. This energy depen- 
dence introduces in turn energy dependence of the moments, as clearly seen 
from equations d6j), even if function F is energy independent. Unfortunately 




(7) 
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this dependence alone would contradict the data since for constant F n mul- 
tiplicity moments decrease with energy (for growing (n)). 

Therefore the source function F has to depend on energy and its mo- 
ments have to win over the decrease generated by the multiplicity growth 
through Eqs.([6]). In Ref.[B] the method of recovering F from the data has 
been discussed, without, however, reference to the recent measurements at 
the LHC. In DPM or QGSM violation of the KNO scaling proceeds by an 
increase of the number of sources (chains, pomerons) with increasing energy. 

Here, rather than constructing a microscopic model of multiparticle pro- 
duction, we choose the explicit form of F(t) and check whether we are able 
to describe multiplicity moments measured by Alice and CMS. To this end 
we choose for F negative binomial distribution (NBD) |10j : 

F(t, k) = ^t k - l e- kt (9) 

which is known to describe relatively well the data at lower energies pQ. 
Distribution (|9p depends on one parameter k, which - as explained above 

- has to depend on W. It is known from the analysis of lower energy data 
that, depending on energy, k ~ 4 - 2 and decreases with increasing energy. 
Let us remind that for k — 1 the probability distribution Pnbd is given 
by geometrical distribution P(n) = (n) n /(l + (n)) n+1 . With increasing k 
(1/k — > 0), the distribution Pnbd is getting narrower tending to the Poisson 
distribution. 

Based on experimental evidence of the wide occurrence of NBD, several 
possible explanations have been proposed in the literature (for review see 
Ref.p]). The NBD has been mostly interpreted in terms of (partial) stim- 
ulated emissions or cascade processes [13]. More recently NBD has been 
derived from the Color Glass Condensate (CGC) approach giving explicit 
prediction for the energy dependence of parameter k at high energies being 
of the order of the LHC energy range [2j. Here, contrary to the lower energy 
trend, parameter k is expected to grow with energy, as it is directly connected 
to the saturation scale which increases with energy. Similar behavior is found 
in String Percolation Model (SPM) [15], where - once percolation is achieved 

- k starts to grow with energy like in the CGC. It is therefore interesting to 
see if the new regime of growing k has been already achieved at the LHC, 
which is one of the motivations behind the present work. 
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For negative binomials 



fc(fc + l)...(fc + m) 
^Wi ~ j^Ti • ^ 



The first equation of ^ gives: 



Using (11) we get for higher moments 

C 2 -l 



C 3 = C 2 (2C 2 - 1) - 



(n> 



C 4 = C- 2 (6^-7<7 2 + 2)-2^~ 4C ' 2 + 1 ■ 6,2-1 



2 ' 



("> (n) 

C 5 = C 2 (24C 2 3 - mCl + 29C 2 - 6) - 2 ™% - U ° l ± 19 ° 2 - 3 

(n) 

UCj - 23C 2 + 9 C 2 - 1 

I \2 / \3 ■ 

(n> (n> 

Let us first observe that for constant C 2 , which is approximately true, at 
least for \r)\ < 0.5 where C 2 ~ 2, higher moments grow with energy. This is 
depicted in Figjl] where long dash (red) line corresponds to constant C 2 in 
two rapidity intervals \rj\ < 0.5 (C 2 ~ 2) and \rj\ < 2.4 (C 2 — 1.54). It is 
clearly seen that for negative binomial distribution used here this growth is, 
however, too slow. For larger rapidity intervals A77 multiplicity is also larger 
dn/drj = 0.755 At] W ' 23 and therefore the inverse powers of multiplicity are 
less important than for smaller rapidity ranges. Moreover for \rj\ < 2.4 second 
moment C 2 ~ 1.54 and therefore the coefficients in front of powers (n)~ m are 
also smaller than for \r)\ < 0.5. Therefore, as seen in FigjTJ for \r)\ < 2.4 in 
the fit where C 2 =const., all other moments are nearly constant as well. 

In order to reproduce the growth seen in the data we have therefore to 
require a mild increase of C 2 with W . Higher moments C m are proportional 
to higher powers of C™ and should therefore grow faster with W with in- 
creasing m. This trend is clearly seen in the data. To this end we choose to 
approximate C 2 by a linear function of log W: 

C 2 = a + b\og(W[GeV}). (13) 
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Figure 1: Multiplicity moments measured by UA5 [TS], Alice [3] and CMS [1] for two 
rapidity intervals in function of energy W = y/s. Long dash (red) lines correspond to 
constant Ci moments equal to 2 and 1.54 for two rapidity intervals |7y| < 0.5 and \r\\ < 2.4 
respectively. Short dash (black) lines correspond to energy dependent as explained in 
the text. 



In order to find parameters a and b we choose to fit C4 rather than C2. 
In both cases \rj\ < 0.5 and < 2.4 moment C4 grows rather fast with 
W having still reasonable errors. Fitting C2 gives usually too slow increase 
of higher moments, whereas fitting C4 reproduces all moments with good 
precision. This is easily seen from Figjl} The parameters of the fit (13) are 
given in Table [l] and the resulting energy dependence of 1/k is plotted in 
Fig. [2j We see that the trend from lower energies continues: k decreases 
with energy but is still rather far from k — 1. For rj < 0.5 and W = 0.9, 7 
and 14 TeV, k = 1.58, 1.25 and 1.18 respectively. Should this dependence 
continue, k — 1 would be reached for W ~ 250 TeV. On the other hand 
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M < vo 


a 


b 




1 


0.5 


1.702 


0.071 




T 


2.4 


0.997 


0.175 




3 


0.5 


1.597 


0.111 




3 


1.0 


1.377 


0.118 




3 


1.3 


1.219 


0.149 



Table 1: Values of parameters a and b of Eq.(13) used in Figs [T] and [3] 



W [TeV] 


10 


14 


10 14 




\v\ 


< 0.5 


\rj\ < 2.4 


(n) 


6.28 


6.78 


31.61 34.15 




1.98 


1.99 


1.70 1.72 


^3 


5.73 


5.81 


4.03 4.18 


c 4 


21.74 


22.30 


12.33 13.11 


c 5 


102.11 


106.02 


46.07 50.39 



Table 2: Predictions for multiplicity moments. 



dependence of k on rapidity is quite pronounced. 
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Figure 2: Growth of l/k used in Fig [I] 
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With parametrizations <ph and (13) we are able to predict the first mo- 



ments for higher energies at which the LHC will be running in the future. 
The results are displayed in Table [2j 

Alice collaboration published results for the multiplicity moments C2 — C4 
for two energies: 0.9 and 7 TeV and three rapidity intervals |?7| < 0.5, 1 and 
1.3. We repeated the same procedure described above for the Alice data 



fitting C4 with the help of Eq.(13). The resulting parameters are collected 
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Figure 3: Multiplicity moments in different rapidity intervals as measured by Alice [5] 
together with NBD fit with C2 parametrized as in Eq.Q and Table [TJ 



in Table [T] and the moments are plotted in Fig|3} We see good agreement of 
NBD fits for all three rapidity intervals. 

To conclude: we have used the convolution model ([T]) with distribution 
of sources given by negative binomial function ^ to fit multiplicity mo- 
ments measured recently by Alice and CMS collaborations at the LHC We 
have shown that convolution model implies that normalized C m multiplicity 
moments decrease with increasing energy as inverse powers of the average 
multiplicity ( 11|12 ) if the distribution function F(t) is energy independent. 



Such a behavior contradicts data. Assuming NBD for F(t) and logarithmic 



growth (13) of C2 moment, we have been able to reproduce the multiplic- 
ity moments over the wide range of energies for different rapidity intervals. 
The input growth of C2 with energy can be easily translated to a decrease 
of the parameter k of NBD function ([9]). This behavior is consistent with 
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lower energies and does not exhibit the change predicted by CGC [H] and/or 
SPM [15] . We also made predictions for higher energies which will be soon 
accessible at the LHC. Unfortunately we are still lacking a microscopic model 
explaining energy dependence of k parameter of the NBD distribution. 
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